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Abstract 

We study the conditions required to make the 2HDM scalar potential stable up to the Planck scale. 
The lightest CP-even scalar is assumed to have been found at the LHC and the alignment limit is imposed 
in view of the LHC Higgs data. We hnd that ensuring stability up to scales > 10^° GeV necessitates the 
introduction of a soft breaking parameter in the theory. Even then, some interesting correlations between the 
nonstandard masses and the soft breaking parameter need to be satisfied. Consequently, a 2HDM becomes 
completely determined by only two nonstandard parameters, namely, tan/3 and a mass parameter, mo, with 
tan /3 > 3. These observations make a 2HDM, in the stable alignment limit, more predictive than ever. 


1 Introduction 


With the discovery of a Higgs-like particle at the Large Hadron Collider (LHC) [1,2], all the parameters of 
the Standard Model (SM) are now known and the fate of the SM is sealed. Taking all the uncertainties of 
the current experimental data into account, it has been found that the SM scalar potential becomes unstable 
somewhere within 10^-10^^ GeV [3-6]. In fact, it has been shown that for the SM scalar potential to be stable 
all the way up to the Planck scale (Mp = 10^^ GeV) the Higgs mass {mh) needs to be in the following range [5]: 

ruh > 129.6 GeV + 2.0 {Mt[GeV] - 173.34 GeV) - 0.5 GeV ± 0-3 GeV, (1) 

where, Mt denotes the top quark pole mass and as{Mz) is the strong coupling constant at the Z-boson mass 
scale. Hence, an SM Higgs boson with mass in the range 124-126 GeV certainly disfavors the possibility of 
having an absolutely stable vacuum up to Mp. As a way out of this vexing situation, it has been suggested 
that while absolute stability of the SM potential might be a tall ask, a metastable vacuum is entirely consistent 
with the current experimental value of the Higgs mass [4,7]. Nevertheless, the problem of vacuum stability in 
the SM remains one of the most discussed topics after the Higgs discovery and often has been taken as a hint 
for the intervention of some new physics. 

Our objective in this paper is to investigate whether the two-Higgs-doublet models (2HDMs) [8] can do better in 
this respect. 2HDMs extend the scalar potential of the SM by adding one extra scalar doublet and therefore, rank 
amongst the simplest of beyond the Standard Model (BSM) constructions. For decades, 2HDMs have attracted 
a lot of attention because the minimal supersymmetry relies on a 2HDM scalar structure. Another attractive 
feature of 2HDMs is that the value of the electroweak p-parameter remains unity at the tree level. But one 
ominous consequence of introducing one additional scalar doublet is that now we will have two Yukawa matrices 
for each type of fermion and diagonalization of the fermion mass matrix will not guarantee the simultaneous 
diagonalization of the Yukawa matrices. As a result, there will be flavor changing neutral currents(FGNG), 
at the tree level, mediated by neutral scalars. This problem was addressed by Glashow and Weinberg [9] and 
independently by Paschos [10]. According to Glashow-Weinberg-Paschos theorem, tree level FCNC can be 
avoided altogether if suitable arrangements are made such that fermions of a particular charge receive their 
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masses from a particular scalar doublet. Usually, a Z 2 symmetry, under which ^ and 02 ^ “ 02 , is 
employed to achieve this. Proper assignments of the charge to different fermions then complete the purpose. 
Labeling as 02 the doublet that couples to the up-type quarks, the following four conventional variants of 
2HDMs emerge from the Z^ charge assignments to the fermions: 

• Type I: All quarks and leptons couple only to the doublet 02 ; 

• Type II: 02 couples to the up-type quarks and 0i couples to down-type quarks and charged leptons; 

• Type X or lepton specific: All quarks couple to 02 , while 0i couples to the charged leptons; 

• Type Y or flipped: Up-type quarks and charged leptons couple to 02 and all down-type quarks couple to 

01 . 


Currently the measured values of Higgs signal strengths into different decay channels are consistent with the 
corresponding SM expectations. In anticipation that the data will continue to agree with the SM with increasing 
accuracy in the upcoming experiments, those BSM scenarios which can deliver an SM-like Higgs in some 
limiting case, will hold the upper hand in future survival race. It is quite well known that 2HDMs promise an 
alignment limit [11-13] when an SM-like Higgs can be recovered in the form of the lightest CP-even scalar in the 
2 HDM particle spectrum. Although vacuum stability constraints in 2HDMs have been studied previously both 
before [14] and after [15] the Higgs discovery, the distinct implications of the alignment limit in this context 
have never been emphasized before. Motivated by the LHC Higgs data, in this paper we concentrate exclusively 
on the alignment limit and explore the consequences of demanding the stability of the 2HDM potential all the 
way up to the GUT and Planck scales. In the process we have uncovered many new and interesting features. 
We have found that the requirement of stability of the 2HDM scalar potential compels us to introduce a soft 
breaking parameter and at the same time, entails a strong correlation between the soft breaking parameter and 
the other nonstandard masses. Accordingly, in this ‘stable alignment limit’ a 2 HDM is completely determined 
by only two nonstandard parameters: (i) tan 0 (= U 2 /U 1 ) which is the ratio of the two vacuum expectation 
values (vevs) and, (ii) a mass parameter (mo). How the requirement of high scale stability of the 2HDM scalar 
potential in the alignment limit leads us to this intriguing conclusion constitutes the central theme of our paper. 

The paper is organized in the following way: in Section 2 we describe the model and the various constraints 
that we use in our study. We present our results in Section 3. In Section 4 we summarize important findings 
and draw our conclusions. 


2 The scalar potential 


The general 2HDM potential with Z 2 symmetry under which 0i ^ 
V = + TO22</>2</>2 - (ml24>\4>2 + h.C.^ + 


01 and 02 ^ — 02 , is usually written as 



+A3 



+ A4 



( 2 ) 


where the term proportional to m ^2 breaks the Z 2 symmetry softly. For simplicity, we assume that all the 
parameters of the potential are real so that CP symmetry is not explicitly broken in the scalar sector. To 
minimize the potential, we express the doublets as 


and define the minimum as follows: 
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Note that, although the potential of Eq. ( 2 ) is explicitly CP conserving, the vevs can, in general, be complex 
[16,17]. Now the minimization conditions will read as follows: 
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2 {mj^xi - mj2X2) + Mxi{xj + yf) + xi{X 3 + \i){xl + y|) 
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2{ml2X2 - ml^xi) + \2X2{xl + yl) + X2{\:i + Xi){xl + yl) 
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+X5{y2yl + 2xiX2yi - xly2) = 0 . 


(5a) 

(5b) 

(5c) 

(5d) 


It can be easily checked that the choice yi = y 2 = 0 which implies the vevs are real, satisfies Eqs. (5c) and (5d) 
trivially. In this case, Eqs. (5a) and (5b) will take the following simpler forms: 

2 (^ 11 ^! - mi 2 '^ 2 ) +'^1 + (As + A 4 + A 5 )n|} = 0 , ( 6 a) 

2 (m 22 '^ 2+'^2 {A 2 n| + (As + A 4 + A 5 )ni} = 0. ( 6 b) 

Thus the assumption of real vi and V 2 is consistent with the minimization conditions of Eq. (5) and in what 
follows, we shall work under this assumption. 
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When both the doublets receive vevs, the Z 2 symmetry is broken spontaneously too and we can rewrite the 
2HDM potential in the following form [18]: 


V = 


A ^0101 -+ P2 ^ 02^2 -+ ^3 ^0101 + 02^2 — 

+/^4 - ((/>l 02 )( 4 ^l)} +/^5 (Re 


t / viV2y 

A-—j 
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In comparison with the previous case, in this notation plays the role of the soft breaking parameter. Eor 
easy understanding of the future results, we will often switch between the parametrizations of Eqs. ( 2 ) and (7). 
Therefore, a comparison between the parameters of Eqs. ( 2 ) and (7) will be in order. It should be emphasized 
that unlike Eq. (2), Eq. (7) manifestly assumes that both the doublets acquire vevs so that mh and 17122 
( 2 ) can be traded for vi and V 2 (using Eqs. ( 6 a) and ( 6 b)) followed by suitable rearrangements of the quartic 
parameters to obtain Eq. (7). Since, in this paper, we restrict ourselves only to the non-inert cases where tan/3 
is nonzero and finite, above two parametrizations are equivalent to us. The connections between the two sets 
of parameters are given by the following relations: 


^11 = -{Pi4 ^ ; 

= -{P 2 VI ^ ; 

2 /^5 

mi 2 = yUiU 2 ; 


A 4 


Pb + Pq 

2 


- ^4 ; 


Ai — 2{pi + ^ 3 ) ; 

A2 = 2(^2 + ^ 3 ) ; 
A 3 = ( 2^3 + Pa) ; 
7 . p5 - Pe 


( 8 ) 


At this point, it is interesting to note that in the parametrization of Eq. ( 2 ), the combination mi 2 /(sin;d cos P) 
instead of m ^2 hself controls the nonstandard masses [11]. Thus the relation between m ^2 /^5 Rd- ( 8 ) 
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suggests that is a better parameter for tracking the effect of soft breaking. In passing, we note that in the 
limit A 5 = 0 (in Eq. (2)) or equivalently (in Eq. (7)), the symmetry of the 2HDM potential is enhanced 

from softly broken Z 2 to softly broken U{1) under which (j)i and <p 2 This U{1) symmetry will 

be relevant in our future discussions. 


To get the mass eigenstates, we expand the scalar doublets around their vevs as follows: 
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Our assumption that CP is a good symmetry of the scalar potential actually allows us to define neutral scalar 
eigenstates that are also eigenstates of CP. In total, there are five physical eigenstates: two CP-even scalars 
(h, i7), one CP-odd scalar (^4) and a pair of charged scalars {H^) along with three Goldstones (G^, G^) which 
will be absorbed to give masses to the SM gauge bosons (IT^,Z). The rotations that lead us to the mass 
eigenstates are given below: 
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where, c^(a) = cos/3(a) and = sin/3(a). The mixing angle of the CP-even sector is defined through the 
following relation: 


tan 2 a = -- 


2(/33 + ^)VlV2 

/?2C| + (/?3 + f)(v! 



( 11 ) 


It is now instructive to count the number of free parameters in the scalar potential. Note that, Eqs. ( 2 ) and 
(7) both contain eight free parameters. In the notation of Eq. (7), these are ui, V 2 and six Pi couplings. We 
can trade vi and V 2 for v = v\ ^ v\ and tan/3. Except for /^s, all other /3 parameters may be traded for 
four physical scalar masses (mh^mH^TriA and mfj+) and the angle, a. The equivalence of these two sets of 
parameters is demonstrated by the following relations: 
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( 12 a) 

( 12 b) 

( 12 c) 

( 12 d) 

( 12 e) 


Among the eight redefined parameters that appear on the RHS of Eq. ( 12 ), not all are unknown. We already 
know V = 246 GeV and under the assumption that the lightest GP-even Higgs is what has been found at the 
LHC, ruh ~ 125 GeV is also known. 


As a next level of simplification, we recall that the experimental values of the Higgs signal strengths into different 
decay modes are showing increasing affinity towards the SM expectations [19,20]. This encourages us to work 
in the alignment limit [ 11 ]: 

/3 - a = I > (13) 
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which means, h will have the exact SM-like tree-level couplings with the fermions and the vector bosons. The 
recent global fits of the LHC data in the 2 HDM context [21-27] certifies that Eq. (13) is indeed a reasonable 
assumption. In what follows, we will work exclusively in the alignment limit. Therefore, in this limit, we are left 
with five unknown parameters tyia^ ^5 arid fail P) which will be constrained from the requirement of 

high scale stability of the 2HDM potential. We will see that this requirement in association with the constraints 
from electroweak T-parameter entail strong correlations between most of the remaining parameters making the 
2HDM more constrained than ever. 


2.1 Theoretical constraints from vacuum stability and unitarity 


First we have to ensure that there is not any direction in the field space along which the potential becomes 
infinitely negative. In the parametrization of Eq. (2), the necessary and sufficient conditions for the potential 
to be bounded from below read [28-31] 


VSCl ; 

V 

o 

(14a) 

VSC 2 : 

A 2 > 0 , 

(14b) 

VSC3 ; 

A3 + V A 1 A 2 > 0, 

(14c) 

VSC4 ; 

A 3 + A 4 — A 5 + a/ A 1 A 2 > 0. 

(14d) 


The ^'-matrix eigenvalues that will be constrained from unitarity of the scattering amplitudes are also listed 
below [32-35]: 


af 


(15a) 

af 

= + ^ 2 ) ^- MY + , 

(15b) 

«3 

= + M) ^{M - MY + , 

(15c) 


= A 3 + 2 A 4 ± 3 A 5 , 

(15d) 

bt 

= A 3 ± A 5 , 

(15e) 

bf 

= A 3 ± A 4 . 

(15f) 


The requirement of tree unitarity then restricts the above eigenvalues as follows: 

\af\,\bf\<16Tr. (16) 


We impose that the inequalities (14) and (16) should be satisfied at all energies between the electroweak 
and Planck scales. The renormalization group (RG) equations that we use to calculate the lambdas at any 
intermediate energy scale are given in Appendix A. 

Additionally it is also important to check whether the minimum defined by vi and V 2 is a global minimum or 
not. The condition for the global minimum, in the notation of Eq. (2), is given by [36,37], 


^ = mi2 




> 0 . 


(17) 


2.2 Experimental constraints 

In addition to the theoretical constraints mentioned above, we also take the following experimental facts into 
account. 
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• There is a very strong lower limit on mfj+^ for Type II models, arising mainly due to the excellent 
agreement of the b ^ sj branching ratio with the SM prediction. Because of this we take 171^+ > 300 
GeV [38,39] for Type II models. However, if tan/3 > 1, there is no such bound for Type I models from 
flavor data [38]. Therefore, for Type I models, we only consider the direct search limit mH+ > 80 GeV [40]. 
At this stage, it should be noted that since the Yukawa couplings in the quark sector are the same for 
Type II and Y models and the top-Yukawa gives the dominant fermionic contribution in the RG equations, 
whatever we comment on Type II models in this paper will eventually hold for Type Y models too. The 
same is true for Type I and X models. 

• The oblique T-parameter can restrict the splitting between the heavy scalar masses. In the 2HDM 
alignment limit, the expression for the new physics contribution to the T-parameter can be expressed 
as [41,42] 


AT = 


with, 


iDTTSin 

[ 0 ior X = y . 


Taking the new physics contribution to the T-parameter as [43] 

AT = 0.05 ±0.12, 


(18) 

(19) 


( 20 ) 


we use the 2(7 uncertainty range around the mean value for our numerical constraints. Note that the 
function F{x,y) is symmetric under x ^ y and is sensitive only to the difference \x — y\. Thus AT = 0 
when either m^+ = tyih or m^+ = In these cases AT puts no constraints on the heavy scalar masses. 
But if, for some reason, rriH ~ ttia then AT severely restricts the splitting between charged and neutral 
scalar masses. 


3 Numerical analysis and results 


For the purpose of analysis, we choose the set {t;, tan a, m/^,, mn^ tua, which appears on the RHS of 

Eq. (12) to be our independent set of eight parameters to describe the 2HDM scalar potential. Among them, 
we set V = 246 GeV and a = P — n{2. In the case of an exact Z 2 symmetry we also set yds = 0 and perform a 
random scan for the rest of the parameters in the following ranges: 

tan/3 G [0.1,10] , rrih G [124,126] GeV , rriH G [rrih, 1500] GeV , rriA G [0,1500] GeV , 

mH+ G [80,1500] GeV (for Type I models), 1 x 1 ^+ G [300,1500] GeV (for Type II models). (21) 

When the Z 2 symmetry is softly broken, we find it convenient to introduce the mass parameter 


mo = 



( 22 ) 


and vary mo in the range [0, 1500] GeV. In the case of a softly broken 1/(1) symmetry, we impose mo = m^- 

Next we use Eq. (12) to convert our set of parameters into the set of ^iS followed by Eq. (8) to convert the /3^s 
again into A^s. We then use the RG equations for A^s given in Appendix A to compute them at any intermediate 
energy scale and check that the unitarity and stability conditions hold all the way up to the Planck scale. We 
set the top quark pole mass at 173.3 GeV [44] for our numerical analysis. Our observations for different cases 
appear in the following subsections. 
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Figure 1: (Exact Z 2 ) Allowed points in the Log]^o(^^^ ci'^d mA-TnH+ planes for the 2HDM potential 

to he stable up to Auv- The points in different colors correspond to different choices for Auv which appear 
in the legends. 


3.1 Exact Z 2 

In this case setting = 0, we scan the rest of the parameters in the range specified previously. The result has 
been displayed in Fig. 1 . A couple of noteworthy features that emerge from this figure are given below: 

(i) The value of tan/3 is bounded and the bound depends on the energy scale, Auv, up to which stability is 
demanded. As we increase Auv, the allowed parameter space shrinks continuously. It should be noted 
that a limit on tan (3 for exact Z 2 symmetry is not surprising at all and is a direct consequence of unitarity 
and stability at the electroweak scale [45,46]. 

{ii) It is found that, for the non-inert cases, neither Type I nor Type II models can be absolutely stable all the 
way up to the Planck scale when the Z 2 symmetry is exact in the potential. This result is in agreement 
with the previous analysis of Ref. [15] in the context of exact Z 2 symmetry. We have checked for each order 
of magnitude in Auv and noticed that Type I models can remain stable up to a maximum of 10^ GeV 
whereas Type II models can be stable only up to 10^ GeV. This difference between Type I and Type II 
models may be understood by noting that stability up to 10^ GeV for Type I models requires a light 
charged scalar, niHr ^ 180 GeV (see lower left panel of Fig. 1 ), which is not allowed for Type II models 
from b ^ sj. 

Now that we know 2HDMs with exact Z 2 symmetry fail to maintain stability up to the Planck scale, it is time 
to investigate whether the introduction of a soft symmetry breaking parameter can improve the situation. 


7 











































1400 
1200 
^1000 
O 800 
^ 600 
400 
200 



Figure 2: (Softly broken Z 2 ) Allowed points in tan^^-mo plane for the 2HDM potential to he stable up to 
Auv* The points in different eolors eorrespond to different ehoiees for Auv whieh appear in the legends. 


3.2 Softly broken Z 2 

Here we study the effects of nonzero For convenience we have encoded the information of into the mass 
parameter, mo, through Eq. (22). The same analysis as in the case of exact Z 2 symmetry is then performed for 
three representative choices of Auv? namely, 10^^, 10^^ and 10^^ GeV. We exhibit our results in Figs. 2 and 3. 
Some intriguing features that emerge from these figures are listed below. 

(а) From Figs. 2 and 3, we note that it is indeed possible for a 2HDM potential to remain stable all the way 
up to the Planck scale in the presence of a soft symmetry breaking parameter. But certain conditions need 
to be satisfied for this to happen. These conditions which we describe below, together define the stable 
alignment limit for 2HDMs. 

( б ) From Fig. 2 we see that for the 2HDM potential to be stable up to Auv? there must exist a lower bound 
on tan/3. For example, when Auv = 10^^ GeV this bound reads tan^d > 3. We also note from Fig. 2 that 
there is a lower limit on mo (or equivalently too. For Type I models we read mo ^120 GeV whereas for 
Type II models we obtain mo ^ 280 GeV. For the energy scales that we are concerned with, we observe that 
these limits on mo for Type I and II models are essentially independent of Auv- We emphasize that this 
is one of the important differences between our result and that of Ref. [15] where the stability of 2HDMs 
with softly broken Z 2 symmetry was analyzed including tan (3 = 2 without realizing that such a low value 
of tan/3 is forbidden when the alignment limit is exact. 

(c) From Fig. 3 we notice that the requirement of absolute stability up to Auv entails a strong correlation 
between mo and the other nonstandard masses. In fact, a closer scrutiny of the plots reveals that all the 
nonstandard scalar masses have to be nearly degenerate with mo- 

{d) Observations ( 6 ) and (c) together lead us to the conclusion that, in the stable alignment limif a 2 HDM 
becomes completely determined by only two nonstandard parameters, namely, tan (3 and a mass parameter 
(mo). Additionally, both of them are bounded from below: 

tan/3 > 3 , mo > 120 GeV (Type I), mo > 280 GeV (Type II). (23) 

(e) We have also checked that above conclusions do not crucially depend on the input parameters, especially 
the top-quark mass. To be quantitative, instead of the central value of 173.3 GeV, if we consider the 2 cr 
lower limit of the top-quark pole mass, 171.8 GeV [44], the lower limit on tan/3 changes to tan^d > 2 . 8 . 

We feel that an a posteriori explanation of most of the above features is possible, at least on a qualitative 
level. In the following we present, one by one, the steps of our argument that helps us apprehend different 
characteristics of Figs. 2 and 3. 
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Figure 3: (Softly broken Z 2 ) Allowed points in niQ-niH, mo-mH+, mo-iriA planes for the 2HDM potential 
to he stable up to Auv- The points in different eolors eorrespond to different ehoiees for Auv whieh appear 
in the legends. 


( 1 ) First we note from Eq. (41) in association with Eqs. (42e), (43e), (44e) and (45e) that the evolution of 
A 5 is proportional to itself. This is not surprising because, as mentioned earlier, in the absence of A 5 
the symmetry of the scalar potential is enhanced to a global U{1). Thus if we start with A 5 = 0 at the 
electroweak scale, it will remain zero at all energy scales. But any initial nonzero value of A 5 will cause it 
to grow with energy. This may eventually jeopardize the vacuum stability condition (14d) at high energy 
scales. Therefore, it will be no wonder if, by demanding high scale stability of the 2HDM potential, we are 
led to the U{1) limit, A 5 ~ 0. Evidently, this limit will have the following implication on the masses (using 
Eq. (12e) and remembering that ^ 0 implies /3^ ^ /3q): 

( 24 ) 

which has been clearly depicted by the first row of plots in Eig. 3. 
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( 2 ) Next important thing to note is that unitarity and stability conditions at the electroweak scale imply, among 
other things, [46,47] 


0 < {mjj — mo)(tan^ f3 + cot^ f3) + 2 m^ < —-— . (25) 

This inequality has been plotted in Fig. 4 for three different values of tan/3. It is obvious that for tan^d 
away from unity, the inequality (25) renders a degeneracy between ttih and rriQ. This explains the second 
row of Fig. 3. 



Figure 4: Allowed regions in the mo-iriH 
plane from the inequality (25) for three differ¬ 
ent values of tan p. 



Figure 5: Allowed region at 95% C.L. from 
the experimental limit on AT under the as¬ 
sumption mn = 'tua = ^ 0 - 


(3) Collecting ( 1 ) and ( 2 ) together we get ra\j ~ m\ ~ m^. Feeding this information to Eq. (19), we obtain 
the following expression for AT: 


AT = 


1 


87r sin^ 0 . 




F{ 


m 


H+^ 


m\ 


D- 


(26) 


Since F (m^+,mo) restricts the splitting — mgl, the experimental limit on AT will impart the de¬ 

generacy between m^+ and mo as shown in Fig. 5. This explains the third row of Fig. 3. 

Thus, up to this point, we have 


^ m%+ , (27) 

which summarizes the essential features of Fig. 3 that the 2HDM can be described by a single nonstandard 
mass parameter. Eq. (27) also explains the lower bounds on mo in Eig. 2 . Since m^+ > 300 GeV for 
Type II 2HDMs, the degeneracy of Eq. (27) destroys the possibility of having a light pseudoscalar for 
Type II models. This result is apparent from the upper right panel of Eig. 3. 

(4) To understand the lower limit on tan (3 we turn our attention to the RG evolution equation of A 2 which is 
given by Eq. (41) in association with Eqs. (42b), (43b), (44b) and (45b). Since we have assumed 02 to be 
the doublet that gives masses to the up-type quarks, A 2 will face the negative pull of the top-Yukawa {ht). 
Note that in the limit of exact degeneracy, m^ = m\ = m|^ = m^+, we have at the electroweak scale, 
(using Eq. (12) followed by Eq. ( 8 ) in the alignment limit) 

2 

Ai = A 2 = A 3 = —^ , and, A 4 = A 5 = 0 . (28) 
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Figure 6: (Softly broken Z 2 ) Running of sta¬ 
bility conditions of Eq. (I 4 ) in the exactly de¬ 
generate scenario assuming nth = 126 GeV. In 
this scenario, these runnings are independent 
of mo. 



Figure 7: (Softly broken Z 2 ) Allowed points in 
plane for Auv = 10^^ GeV in a Type I 
2HDM with mn = mA = mo, and m^r = 
mo + 8. 


Using these, we may simplify the RG equation for A 2 , at the electroweak scale, as follows: 

12 X 2 = I 6 A 2 — 3 ^3^^ g' ^ ^ ^3^^ g' + ‘^ 9 ^ 9 ' ^ + 12h^A2 — 12h^ . (29) 

Notice the striking similarity between the above equation and the SM running of A that appears in Eq. (35). 
One difference is that PA 2 receives additional contributions on the positive side (compare I 6 A 2 with 12 A^ 
in the SM case) due to the presence of extra quartic couplings. On the contrary, ht in the 2 HDM case is 
r\j \/2mtl{v siu/d) which is larger than the SM value of ht ^ \/2mtlv. Therefore, compared to the SM case 
the negative drag of ht is enhanced in 2HDMs. This effect drives A 2 to negative values violating vacuum 
stability condition (14b) at high energies unless sin [3 is large enough to dilute the effect of the term, — 12 h|, 
sufficiently. This is the origin of the lower limit on tan^d. Choosing tan^d = 8.5 for illustration, we have 
displayed this effect in Fig. 6 . There we see how the evolution of the dotted (blue) line marked as VSC 2 
(which is nothing but A 2 ) marginally survives becoming negative at an intermediate energy scale. Take 
tan/3 too low and this line goes below zero spoiling the stability condition, A 2 > 0. It is worth noting that 
when all the nonstandard masses are exactly degenerate, evolutions of the A^s depend only on tan^d and 
not on the mass parameter, mo. 

(5) With exact degeneracy, the lower bound on tan P turns out to be tan > 8 for Auv = 10^^ GeV. But from 
Fig. 2 we see that even lower values of tan/3 are allowed. To understand this, we need to investigate the 
evolution of A 2 in greater detail. Instead of exact degeneracy, we now consider the following limit: 

mo = mA = mH , and, + A . (30) 

In this limit, 

777^ 2 A 2 A 

Ai = A 2 = , A 3 = A 2 H-^ , A 4 =-^ , A 5 = 0 . (31) 

yZ yZ yZ 

Using these, we can rearrange the terms that appear on the RHS in the RG equation for A 2 , to obtain 

11 X 2 = I 4 A 2 + 2 ^A2 + ~ ^ 9 ' ) ^2 + ^ (^^9^ 3 - 9 ' + 2^^^' ^ + 12h^X2 — I2h^ . (32) 

Comparing this with Eq. (29) we see that a positive value of A aids in the positive terms allowing lower values 
of tan /3. But the fact that only very small values of A can be permitted from unitarity [46] and T-parameter 


II 













Figure 8: Chain of arguments that helps us make sense of the main features of Figs. 2 and 3. 

puts a lower limit on tan P anyway. In Fig. 7 we have illustrated how the difference 5 = —tuq ~ A/(2mo) 
depends on tan/3. There we can see that to allow tan/3 < 8 we will need (3 > 0, z.e., m^+ > mo- 


We have summarized our arguments above compactly in the form of a ffowchart in Fig. 8 . Before moving on, 
we intend to make some remarks on the structure of the potential in the stable alignment limit. We note that 
in the limit, = m\ = the potential of Eq. (2) takes the following simple form:^ 

V = +m22</>2</>2 - ^12 + h.c.^ + y + 402 ) • (33) 

This potential now contains four parameters among which two, disguised as v and mh^ have been measured. 
Two other parameters, in the form of tan P and mg, remain to be determined to fix the model completely. Note 
that the symmetry in the quartic terms of Eq. (33) is more than a mere U (1), in fact, the symmetry is enhanced 
to a global U{2) under which (0i, ^ 2 )^ ^ b^(2)(0i, 02 )^- But this U{2) symmetry is explicitly broken in 
the Yukawa terms. Due to this, the structure of the potential in Eq. (33) is not stable under RG. Neither the 
correlation Ai = A 2 = A 3 nor the equality A 4 = 0 is maintained at higher energies. What remains valid at any 
scale is the equality A 5 = 0 because the 1 /( 1 ) symmetry that prevails in the quartic part of the scalar potential 
in the absence of A 5 is also preserved in the Yukawa sector by construction. Therefore, we conclude that when 
high scale stability of the 2 HDM potential needs to be ensured, softly broken U{1) symmetry should be a more 
natural choice, to tackle tree-level ECNC, than the conventional Z 2 symmetry. 


Now we want to check whether the condition (17) for the global minimum holds in the stable alignment limit 
or not. Eirst we remember that, in this limit, tuq is positive (see Eig. 2 ) which implies m ^2 0- We have also 

found Ai ~ A 2 ~ A 3 with tan/3 > 3. These imply 



^Similar potential can also be motivated from an SO(5) symmetry [48]. 


(34a) 
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The remaining factor in Eq. (17) can be evaluated, in the stable alignment limit, as follows: 


m?i 



«^ii - = P2vl - Pivl 

m\ (sin^ (3 — cos^ (3) 
m\ cos^ ;d(tan^ — 1) > 0 


[Using Eq. (8)] 

[Using Eq. (12) with mo ~ mA] 
[•.• tan/3 > 3] 


(34b) 

(34c) 

(34d) 


Therefore we conclude that, in the stable alignment limit, condition (17) is satisfied automatically and hence, 
existence of the global minimum is guaranteed. 


3.3 Softly broken U{1) 

Here we will have ttia = from the beginning. Evidently, when the U{1) symmetry is exact in the potential, 
z.e., mo = 0, H will become the Goldstone boson. We present the result of our analysis, for Type I models, 
in Eig. 9. Although these plots do not provide any new insights, we include them just for completeness. The 
interpretation of Eig. 9 will be similar to that of the corresponding plots in the case of softly broken Z 2 symmetry. 
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Figure 9: (Softly broken U{1)) Allowed points in tan/3-mA, tua-tuh, ruA-'^Hr planes for the 2HDM 
potential to he stable up to Auv- Three different eolors eorrespond to three different ehoiees of Ajjy whieh 
appear in the legends. 


4 Summary and conclusions 


In this paper we have studied the consequences of demanding absolute stability of the 2HDM potential all 
the way up to the Planck scale. In view of the fact that the LHC Higgs data are in conformity with the SM 
expectations, we decided to work in the alignment limit where the lightest CP-even scalar {h) possess SM- 
like tree-level couplings with the standard particles. Although we have explicitly demonstrated our results for 
Type I and H models only, they are equally applicable to Type X and Y models also. We have found that, in the 
alignment limit, requirement of high scale stability puts some remarkable restrictions on the 2HDM parameter 
space. This set of restrictions defines the stable alignment limit for 2HDMs. Some of our most important 
observations are summarized below: 

• 2HDM scalar potential with an exact Z 2 symmetry is unable to maintain stability after 10^ CeV (10^ CeV) 
in the Type I (H) case. To ensure stability up to the Planck scale, Z 2 symmetry needs to be broken softly. 

• By demanding high scale stability in the presence of a soft breaking, we are led to a situation where the 
symmetry of the potential is enhanced from softly broken Z 2 to softly broken U{1). 
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• To have stability up to very high energies (> 10^^ GeV), all the nonstandard masses need to be nearly 
degenerate: mo ~ tua ~ ttih ~ mfj+. Thus, there is only one nonstandard mass parameter that governs 
the 2HDM in the stable alignment limit. 

• The value of tan (3 is bounded from below. 

Therefore, in the stable alignment limits a 2HDM can be fully described by only two new parameters, namely, 
tan /3 and mo with tan f3 >3 and mo ^ 120 GeV (280 GeV) for Type I (II) models. 

Finally, we end the paper on a pessimistic note. One should remember that, in the alignment limit, even 
the self couplings of h are identical to the SM at the tree-level. This means, in this limit, any measurement 
involving the standard particles can smell the presence of nonstandard physics only through loop effects. But 
then the nonstandard particles can be made heavy enough to dilute these effects sufficiently. Thus, our only 
hope to detect the presence of nonstandard scalars beyond the reaches of the direct search experiments is via 
some nondecoupling effects in certain loop induced processes like Higgs to diphoton decay [49]. But with the 
degeneracy, m^ ~ this possibility is also eliminated. Therefore, it appears that even if the LHG data 

continues to agree with the SM predictions, a 2HDM in the stable alignment limit which is very difficult, if not 
impossible, to probe experimentally, can serve as an alternative to the SM for many years to come. 


Acknowledgments: DD thanks the Department of Atomic Energy (DAE), India for financial support. 


Appendix 

A One loop RG equations 


Let us first write down the one loop RG equation (RGE) of the SM quartic coupling (A) [50,51] as follows: 

VX = 12X^ - 3 (sy + 5 '^) A + ^ (sy + 5 '^ + + Uh^X - 12h^ , (35) 

where, ht denotes the top-Yukawa which dominates over the other Yukawa couplings. Eor convenience, we have 
introduce the shorthand V = IGtt^ d(in/x) • 

Now we will present the one loop RGEs of all the relevant couplings (gauge, Yukawa and scalar quartic couplings) 
of 2HDM [8,52-54]. 


Gauge couplings: RGE for the gauge couplings, 

= 

-75? , 

(36a) 

Vg = 

-35^ 

(36b) 

Vg' = 

7/. 

(36c) 


Yukawa couplings: The initial values of top (ht), bottom (h^) and tau (h^-) Yukawa couplings at electroweak 
scale (Mt) for the different Yukawa structures are given by. 


_ 


vsingi^ 37T^s{Mt)} , 
V sin (3 ’ 


Type I: 


ht{Mt) 

hb,T{Mt) 
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(37a) 





Type II: 


(37b) 


ht{Mt) 

hb,r{Mt) 


V cos /3 ’ 


where, as{Mt) denotes the strong coupling constant at top quark pole mass. The corresponding RGBs for 
Type I Yukawa structure are given by 


Vht 


Vhi) 


Vhr 


ht + -h^ + -hi + , 

hb (^CLd + -hi + -hi + hl^ , 
hj- (aQ-\- 3hl + 3hl ^ 


(38a) 

(38b) 

(38c) 


For Type II Yukawa structure, the RGBs take the following form: 


where. 


Vht 

Vhb 

Vh^ 


ht[cLu^ -hi + -hi 


hb ( ad^ -hi + -hi + hi 


hj- { clq + 3hl —hi 


— 


o^d = 


O 2 ^2 /2 

- 85 . -19 

O 2 ^2 ^ /2 

-8^. -19-^2^ 


= -T^ “ 


15 ,2 


(39a) 

(39b) 

(39c) 

(40a) 

(40b) 

(40c) 


for both Type I and II models. 


Scalar quartic couplings: The RGBs for the five quartic couplings that appear in Bq. (2) are given by 


V\t = px, + G, + i/, , (i = 1, 2,3,4, 5), (41) 

where, P\. and Gi are independent of the Yukawa structure of the model and are as follows: 

Pxi = 12A^ + 4 A 3 + 4 A 3 A 4 + 2 A 4 + 2 A 5 , (42a) 

Px2 — I 2 A 2 + 4 A 3 + 4 A 3 A 4 + 2A| + 2 A 5 , (42b) 

— (^1 + ^ 2 ) (bA 3 + 2 A 4 ) + 4 A 3 + 2 A 4 + 2 A 5 , (42c) 

Pxa — 2 A 4 (Ai-f-A 2 ) + 8 A 3 A 4 + 4A| + 8 A 5 , (42d) 

j3x^ = A 5 (2Ai + 2 A 2 + 8 A 3 + I 2 A 4 ) , (42e) 

and, 

Gi = - ^3^^ ^ g' + ‘^g^g' ^ — 3Ai {3g^ ^ 9 ' ^ ^ (43a) 

G 2 = - ^3^^ ^ g' + ‘^g^g' ^ — 3 A 2 {3g^ ^ 9 ' ^ ^ (43b) 

Gi = \ ( 35 ^ + 9 '^ - 25 ^ 5 '') - 3 A 3 ( 35 " + , (43c) 
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(43d) 

(43e) 


G4 = ^g^g' — 3A4 (^ 3 g^ g' ^ ^ 

G5 = —3A5 ^ 3 g^ + g^ ^ • 

The expressions for i7^s, however, depend on the Yukawa structure of the model. For Type I models, these are 
given by. 


w 

= 0, 

(44a) 

H 2 

= 4 A 2 {3h^ + 3h^ + h^) — (l2h^ + f+ 4h^) , 

(44b) 

w 

= 2A3 {3h^ + 3h^ + h^) , 

(44c) 

w 

= 2A4 {3h^ + 3h^ + h^) , 

(44d) 

w 

= 2A5 {3h^ + 3h^ + h^) . 

(44e) 


For the Type II Yukawa structure, we have 


Hi 

= 4Ai {‘ihl + hi) - {l2hl + 4hl) , 

(45a) 

H 2 

= I2A2/1I - 12 hl, 

(45b) 

W 

= 2 X 3 {3hl +3hl + hi) - 12 hlhl , 

(45c) 

W 

= 2A4 (3/1^ 3hl hi) ^2^h^hl , 

(45d) 

w 

= 2A5 (3/i| + 3/1^ + hi) . 

(45e) 
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